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ABSTRACT 

The intermediate polar HT Cam is unusual in that it shows no evidence for dense 
absorption in its spectrum. We analyse an XMM-Newton observation of this star, 
which confirms the absence of absorption and shows that the X-ray spin-pulse is 
energy independent. The modulation arises solely from occultation effects and can be 
reproduced by a simple geometrical model in which the lower accretion footprint is 
fainter than the upper one. 

We suggest that the lack of opacity in the accretion columns of HT Cam, and 
also of EX Hya and V1025 Cen, results from a low accretion rate owing to their being 
below the cataclysmic variable period gap. 

Key words: accretion, accretion discs - stars: individual: HT Cam 
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1 INTRODUCTION 

A characteristic feature of intermediate polars (IPs) - in¬ 
teracting white-dwarf/red-dwarf binaries with a magnetic 
primary - is that their X-ray emission shows a clear mod¬ 
ulation at the white-dwarf spin period (see Patterson 1994 
for an overview of these objects). 

Understanding these pulsations, however, is not straight 
forward since many factors can contribute to the modula¬ 
tion. These include opacity in the X-ray emitting accretion 
regions, and opacity owing to infalling material if it passes 
through the line of sight. Further, they include occultation 
as the accretion regions pass over the white-dwarf limb, in¬ 
cluding the effect of possible asymmetries between the two 
accreting poles. Thus IP spin pulses can show complex pro¬ 
files, as exemplified by FO Aqr (e.g. Beardmore et al. 1998) 
and PQ Gem (e.g. Mason 1997). 

It is thus useful to study systems where we can under¬ 
stand the pulse profile in terms of simple geometry alone. 
One such IP is XY Ari, where an eclipse provides geomet¬ 
rical information (Hellier 1997). Recent work on HT Cam 
suggests that it has an unusually simple spin pulse, since 
there appears to be very little opacity affecting the X-ray 
emission (de Martino et al. 2005), and thus no complica¬ 
tions owing to absorption by accretion curtains. 

In this paper we analyse XMM data on HT Cam and 
produce a geometric model of the pulse, leading to con¬ 
straints on the accretion geometry in this system. 
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2 OBSERVATIONS 

HT Camelopardalis (RX J0757.0-1-6306, hereafter HT Cam) 
was observed by the XMM-Newton satellite (Jansen et al. 
2001) for ~40 ks on 2003 March 24. AAVSO data reveal 
that the system was quiescent at V ~ 17. The EPIC 
MOS (Turner et al. 2001) and PN (Striider et al. 2001) 
instruments were operating in full frame mode, observ¬ 
ing through the Medium filter. Data from these three in¬ 
struments were extracted from circular regions of radius 14 
arc sec, centred on the source. An annulus around this was 
used to estimate the background. The OM camera (Mason 
et al. 2001) was also in operation, collecting 17.5 ks of data 
through each of the B and UVW2 filters. The data were 
analysed using the xmm-SAS software, v6.0.0. 

The lightcurves are shown in Fig. 0 The X-ray data 
show clear modulation at the 515-s white-dwarf spin period. 
Fourier analysis confirms the spin pulse but shows no orbital 
or beat modulation. In the B band there is a ~ 5 per cent 
spin-period modulation, but there is no detectable pulse in 
the UV-band, up to a limit of ~ 14 per cent. The mean X- 
ray count-rate in the pn detector is 1.05 counts s“^, which 
corresponds to a 0.2-12 keV flux of 6.5x 10“^^ erg cm“^ s“^; 
for comparison the 0.1-2.4 keV flux observed by ROSAT was 
1.87x10“^^ erg cm“^ s“^ (Tovmassian et al. 1998). 


3 SPECTROSCOPY 

X-ray emission in an IP comes from multi-temperature plas¬ 
mas beneath the accretion shock. To model the spectrum of 
HT Cam we used mekal components in xspec, absorbed by 
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Figure 1. Lightcurves of HT Cam in the X-ray (EPIC-pn; 0.2-12 keV), B (3900-4900 A) and UV (1800-2250 A) bands. The bins are 
100 s, and the 515-s spin-pulse can be clearly seen in the X-ray band. Typical errors are shown. 


both simple and partial covering absorption. To reduce the 
effects of cross-calibration uncertainties the mekal normal¬ 
isations were permitted to optimise independently for the 
three EPIC instruments. All other parameters were kept the 
same for the three detectors. Using two mekals of different 
temperatures gave a of 1530 (Xy=1.20); adding a third 
MEKAL reduced x^ to 1428 (xS = 1.12), and a fourth gave a 
X^ of 1331 (xS = 1.05). The addition of a fifth mekal did not 
improve the fit. 

The CEMEKL model in xspec reproduces mekal emis¬ 
sion along a continuous range of temperatures, where nor¬ 
malisation varies according to {T/Tmax)°'■ We tried replac¬ 
ing our four discrete mekals with a CEMEKL, but this wors¬ 
ened the fit to a x^ of 1385 (xy=1.08). 

In all of the above fits, the column density of the par¬ 
tial absorber tended to zero. We thus removed it, with no 
affect on the fit quality. We also tried adding a soft black- 
body component, as is seen in some IPs (e.g. de Martino et 
al. 2004) but this gave almost no improvement in £t qual¬ 
ity (Ax^ ^ 4). The EPIC-pn spectrum plus the four-MEKAL 
model are shown in Fig.|^ and the parameters for the model 
are given in Table Q Note that, as the parameters are ex¬ 
pected to vary on the spin cycle, these may be weighted 
averages. 


4 THE SPIN PULSE 

The spin-pulse profile of HT Cam is shown in Fig.|^ We can 
describe the X-ray pulse as a flat maximum covering phases 
^0.7-1.1 and a dip covering ~ 0.2-0.6. To investigate the 
cause of this dip we extracted spectra from the two phase 
regions and htted the four-MEKAL model to both simulta¬ 
neously. Initially we allowed all parameters to vary between 
the two regions, which gave a x^ of 1750 (x2=l-04)- This 
was almost unchanged (x^ = 1753, x2=l-04) if the absorp¬ 
tion was tied between the regions. We tried adding a partial 
absorber to the phase-minimum spectrum, but it made no 
difference to the fit, with the column density tending to zero. 
As a check we tried forcing the emission to remain constant 
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Figure 3. The spin pulse-profile of HT Cam, with typical errors. 
The upper panel shows the mean 0.2-12 keV flux from all three 
EPIC cameras. The softness ratio is defined as (0.2—4)/(6—12) 
keV. The B band is as Fig. 1. Phase zero is defined from the 
ephemeris of Kemp et al. (2002). 


across spin-phase, varying only the absorption. This gave 
a poor fit (x^=1969 x2=l-15). We thus conclude that the 
spin-pulse is not caused by absorption variations, but must 
result from variations in the underlying X-ray emission. 

The softness ratio shows almost no modulation on 
the spin-period, implying that the modulation is energy- 
independent. We thus fixed the absorption and the ratios 
of the emission components across spin phase, attempting 
to reproduce the observed spectral changes by reducing the 
flux during the dip but keeping the ‘colour’ constant. This 
gave a good fit (x^=1780, x^ = l-05) which suggests that 
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Figure 2. The EPIC-pn spectrum of HT Cam, fitted by a four-MEKAL model (see text). Note the prominence of the thermal Iron-K 
emission complex around 6.7 keV. The residuals provide evidence for a 6.4-keV component of cold iron, having an equivalent width of 
70 eV. 


the modulation is indeed energy independent within XMM's 
pass-band. 

In an IP the accretion disc is magnetically truncated, 
with material flowing along the magnetic field lines in large 
‘accretion curtains’. In many systems these curtains peri¬ 
odically pass through our line of sight, causing a promi¬ 
nent absorption dip (e.g. AO Psc, Hellier, Cropper & Mason 
1991; FO Aqr, Beardmore et al. 1998; V1223 Sgr, Taylor et 
al. 1997). The low absorption in our spectroscopy and the 
energy-independent nature of the spin-pulse mean that this 
cannot be happening in HT Cam. Further, they suggest that 
the post-shock accretion columns themselves have an unusu¬ 
ally low opacity, so that we see no difference in the spec¬ 
tra when viewing them at different angles. This contrasts 
with, for example, V405 Aur, where the accretion curtains 
never cross our line of sight, but the changing aspect of 
the columns themselves gives hardness-ratio changes on the 
spin-period (Evans & Hellier 2004). 

We suggest that, in contrast, the spin pulse in HT Cam 
arises solely from simple occultation effects as the accretion 
columns pass over the white-dwarf limb as it rotates. This 
allows us to investigate the X-ray spin pulse with a simple 
geometrical model. 


5 MODELLING THE EFFECTS OF 
OCCULTATION 

We constructed a model which consists of arc-shaped ac¬ 
cretion regions of zero height above the white-dwarf surface 
(since the majority of the X-ray emission in an IP occurs 
close to the white-dwarf surface, see Cropper, Wu & Ram¬ 
say 2000). The shape and location of the arcs was deter¬ 
mined by assuming a disc disrupted at rmag, and by tracing 
field lines from there to the white-dwarf surface, assuming a 
magnetic dipole inclined to the spin axis by an angle 5 and 
offset from the white-dwarf centre in an arbitrary direction. 
A schematic diagram of this model is shown in Fig. 2] Since 
accretion rate will likely vary with azimuth, being greatest 
at the point to which the magnetic pole points, we assume 
a brightness varying as 

B(7) = Acos= (1) 

where 7 is the magnetic longitude of a point on the footprint, 
between end-points ±A 7 , and A is a brightness coefficient, 
which can take different values for the upper and lower re¬ 
gions. By default the upper accretion lay between ±90° and 
the lower one covered 90-270°. The star was then rotated 
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Component 

Parameter 

(Units) 

Value 

Error 

Absn. 

UH (cm“2) 

5.35x10^0 

(±0.28, -0.45) 

Mekal 

kT (keV) 

0.172 

(±0.029, -0.033) 


Abundance 

0.622 

(±0.103, -0.104) 


Norm (Ml) 

4.11x10-® 

(±8.05, -4.11) 


Norm (M2) 

7.63x10“® 

(±9.94, -4.67) 


Norm (PN) 

1.02x10-* 

(±6.24, -3.32) 

Mekal 

kT (keV) 

0.678 

(±0.022, -0.022) 


Norm (Ml) 

1.52x10-* 

(±0.33, -2.67) 


Norm (M2) 

1.46x10“* 

(±0.32, -2.54) 


Norm (PN) 

1.97x10“* 

(±6.53, -0.14) 

Mekal 

kT (keV) 

1.99 

(±0.29, -0.26) 


Norm (Ml) 

5.45x10“* 

(±1.69, -2.09) 


Norm (M2) 

6.88x10“* 

(±0.93, -1.72) 


Norm (PN) 

5.59x10“* 

(±1.38, -1.33) 

Mekal 

kT (keV) 

13.4 

(±2.6, -2.1) 


Norm (Ml) 

2.25x10“® 

(±0.67, -0.16) 


Norm (M2) 

2.00x10“® 

(±0.17, -0.10) 


Norm (PN) 

2.06x10“® 

(±8.23, -0.12) 


Table 1. Model components and parameters fitted to the phase- 
averaged spectrum of HT Cam (see Section n. The errors are 
given to the same power of ten as the values. 


through one revolution in 100 phase bins and a simulated 
pulse-profile was produced by summing the brightness of 
the visible emitting regions. We ignored projection effects, 
treating the X-ray emission as optically thin. 

Some form of asymmetry between the two accretion 
regions is necessary to produce any modulation, otherwise 
the occultation of the upper pole is compensated for by the 
emergence of the lower pole. We explored three sources of 
asymmetry: offsets of the dipole in different directions; re¬ 
ducing the azimuthal extent of one region compared to the 
other; and reducing the brightness of one region compared to 
the other. We explored these asymmetries over a full range 
of system inclination (i) and the dipole inclination (5), in 
steps of 10°. 

A selection of model pulse-profiles is given in Fig. 
and discussed below. Note that if we reverse the sense of 
the asymmetry between the poles, the pulse-profiles are in¬ 
verted. 


5.1 Offsetting the dipole 

We first offset the dipole along the x-axis (perpendicular to 
the plane of the spin and magnetic axes, see Fig. El. Since 
disappearance of one pole is not coincident with appearance 
of the other, there is a modulation. However it is very dif¬ 
ferent from the data (Fig. El top row), so is not considered 
further. 

If the offset is instead in the y-direction (in the plane 
of the spin and magnetic axes), the lower footprint will not 
appear over the limb at all when i and 5 are low. This re¬ 
sults in a deep dip (Fig. El second row). The flux seen at 
minimum in Fig. El results because the upper pole is never 
totally occulted. When i and & are higher the lower foot¬ 
print does appear, but not until after the upper one has 
begun to be occulted, hence a partial dip is seen. Once the 


lower pole does appear it causes the flux to rise, causing a 
‘hump’ during the dip. 

Invoking a 2 :-offset (i.e. along the spin axis) makes the 
magnetic latitude of the footprints bigger at one pole than 
the other. At low values of i and 5 this prevents the lower 
footprint from ever being seen (Fig. El third row). As i and 5 
increase, the lower footprint does appear, but as it is smaller 
than the upper one it appears and disappears more rapidly, 
which gives rise to a ‘hump’. 

None of the above offsets result in a pulse profile resem¬ 
bling that in HT Cam. 

5.2 Reducing the size of the lower pole 

We tried reducing the longitudinal range of the lower arc, 
keeping the brightness coefficient (A) the same, so that the 
lower pole becomes fainter overall. The difference in total 
brightness between the poles gives rise to a dip. The mid¬ 
part of ingress is flat, since appearance of one pole is can¬ 
celling with disappearance of the other. However, the effect 
of the smaller footprint is confined more to the middle of 
the ingress, so a decline is seen at the beginning and end of 
ingress. 

In the above the total brightness of the lower foot¬ 
print is less than that of the upper. We have explored, in¬ 
stead, changing the brightness coefficient such that the total 
brightnesses are the same. This gives a worse match to the 
data, with a ‘hump’ being seen instead of a dip, and is not 
shown in Fig. El 

Again, the asymmetries considered in this section do 
not produce profiles like that seen in HT Cam. 

5.3 Making the lower footprint fainter 

Keeping the dipole symmetrical, but reducing the brightness 
coefficient for the lower footprint, results in a smooth mod¬ 
ulation as the upper footprint is progressively replaced with 
the fainter lower one (Fig. El bottom row). Changing the 
system and dipole inclinations alters the width and profile 
of the dip. If 5 and i are high the entire lower pole is visible 
for some time, producing a flat-bottomed minimum. 

This model gives the only good match to the data. Good 
fits can be found for a system inclination of ^ 40-70°, and 
a dipole inclination of ~ 20-50°. Here the footprints lie be¬ 
tween ±90° (upper) and 90-270° (lower). To reproduce the 
depth of the dip the lower footprint is 65 per cent as bright as 
the upper one. The magnetic latitude of these regions corre¬ 
sponds to a magnetospheric radius of ^ 1.3-3.6 Rwd- Fig. El 
shows the system with the best-fitting parameters (i = 55°, 
5 = 35°, Rmag = 2.6i?wD) and shows that the accretion 
curtains do not intersect our line of sight in this geometry, 
as required by the lack of absorption in the spectra. 


6 DISCUSSION 

We have found there to be little evidence for opacity in the 
X-ray spectrum of HT Cam, suggesting that the spin pulse 
must result solely from geometric effects. We have thus pro¬ 
duced model pulse profiles for a range of possible geometries. 
We find that only one matches the data, namely one in which 
the lower accretion region is fainter than the upper one. By 
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fine-tuning this model to the data we find that the inclina¬ 
tion of HT Cam is ~ 55°, while the magnetic dipole makes 
an angle of ~35° to the spin axis (see Fig. EJ. 

Analysis of an XMM-Newton observation of HT Cam 
has shown the pulse-profile to be energy-independent. By 
creating model pulse-profiles we show that this can be ex¬ 
plained by the upper pole being periodically occulted by 
the white-dwarf, provided that the lower accretion region is 
fainter than the upper one. Our models also suggest that the 
inclination of HT Cam is ~ 55°, while the magnetic dipole 
makes an angle of ~ 35° to the spin axis. 

The energy independent nature of the X-ray spin-pulse 
in HT Cam and the lack of absorption in its spectrum is 
unusual, since most IPs show spectral variations over the 
spin cycle and prominent, phase-varying absorption. One 
part of the explanation is that, in HT Cam, the accretion 
curtains do not cross our line-of-sight to the X-ray emitting 
regions. But this alone is not sufficient to explain the energy- 
independent X-ray modulation: Evans & Hellier (2004) and 
de Martino et al. (2004) report that the same is true in 
V405 Aur, however the pulse profile of that star is energy 
dependent. In that system, the energy dependence is thought 
to arise from opacity within the accretion columns, which 
changes with the viewing angle. 

The fact that we don’t see this in HT Cam suggests 
that the opacity in the accretion column is low. Only two 
other IPs are known in which a dense absorber {uh > 2 x 
10^^ cm“^) is not needed to fit the X-ray spectrum, namely 
EX Hya and V1025 Cen (see Allan, Hellier & Beardmore 
1998 and Hellier, Beardmore & Buckley 1998, respectively). 
Like HT Cam, these system both lie below the period gap, 
where systems are believed to have accretion rates an order 
of magnitude lower than those above the gap (e.g. Warner 
1987). 

Mukai et al. (2003) came to a similar conclusion for 
EX Hya. They analysed Chandra HETG spectra of several 
IPs above the gap and also of EX Hya. They found much 
more line emission in the latter, which they attributed to an 
optically thin accretion column, probably resulting from a 
low specific accretion rate. 

We thus suggest that the X-ray spectra of IPs depend 
on whether they are above or below the cataclysmic variable 
period gap. Those above the period gap have high accre¬ 
tion rates, resulting in optically thick accretion columns and 
prominent phase-varying absorption when accretion curtains 
obscure the line of sight. Those below the gap have much 
lower accretion rates, optically thin accretion columns and 
show much less absorption in their spectra. 
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Figure 5. Simulated X-ray pulse profiles. In the upper rows the dipole is offset by 0.15 i?wD in the x, y- and 2 : directions. The upper 
footprint covers —90° < 7 < 90° and the lower ones covers 90° < 7 < 270° in all but the ‘footprints’ row. Here the lower footprint covers 
120° < 7 < 240°. In the bottom row the lower footprint is 65 per cent as bright as the upper one. 
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Figure 6. The best-fitting model, as described in the text. Left panel: The data and model lightcurve. Right panel: A schematic view 
of the system. The upper and lower accretion curtains are shown as dashed and dotted lines respectively, while the ring shows the inner 
edge of the disc. We can see that the accretion curtains never cross our line of sight to the footprints. 










